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Ph| Abstract 
P 



QJ , Quark-antiquark bound states are examined in the long-range strong- 



coupling limit with the minimal area law of lattice gauge theory assumed 
as input. Matrix element relations are established which in the effective the- 
' ory obtain dynamical equations equivalent to a formulation of the flux-tube 



C3 



model. 



1 Introduction 

The success of QCD as the correct theory of strong interactions has rested primarily 
upon perturbative calculations. The prediction of hadronic properties is inherently 
non-perturbative and has been much slower to develop. Most off lattice progress 
has turned out to be largely irrelevant to hadron dynamics. An exception has been 
the low quark velocity Wilson loop expansion [|TJ . The aim of this paper is to extend 
this method to arbitrary quark velocity. Our principal result is that under natural 
assumptions QCD provides a framework for meson dynamics that is identical with 
the relativistic flux tube (RFT) model |2|, |3| . 

The RFT model assumes that a chromoelectric flux tube stretches in a straight 
line between the quark and the antiquark in a meson. The tube rest frame energy per 
unit length a contributes to rotational energy, momentum, and angular momentum 
in the meson rest frame. Also, as pointed out by Buchmuller a pure chromoelec- 
tric field in the tube rest frame implies no long range spin-spin correlation. In the 
last few years the RFT model has been demonstrated to be numerically calculable 
for any quark mass case || . The RFT model provides intuitive physical pictures for 
relativistic corrections in agreement with QCD and becomes a Nambu-Goto string 
for high rotational states 0. 

The QCD Lagrangian density is a fundamental object in the theory of strongly 
interacting particles. For a meson it is given by the minimally transformed free 
Lagrangian 



Cfree - £ Q CD = £ {% P ~ ~ -F^Ff , 



(i) 

(2) 



from which follow the Euler-Lagrange equations 



(ip- mj)qj 



D»F, 



9 ^QjTjuQj ■ 







(3) 
(4) 
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These are non-linear and amenable only as a perturbation expansion in small cou- 
pling g. Wilson's observation []5| over two decades ago has led to some of the most 
fruitful work in the non-perturbative regime. The simple picture offered in his area 
law for static quarks is compelling: Contributions to the qq propagator fall off expo- 
nentially in the area swept out by lines of gauge invariant chromoelectric flux joining 
the const it uents's world lines; in this way, widely separated paths are suppressed. 
Subsequent lattice simulations verify the same qualitative behavior for dynamical 
quarks. 

Progress in obtaining a consistent description of confinement for realistic off- 
lattice calculations has been largely uncertain. The variety of semi-relativised con- 
finement models based on the above simple idea is striking and is partly due to the 
limited theoretical input available beyond the static limit. Notable exceptions are 
in works of Eichten and Feinberg, Gromes and later those of Brambilla and Prosperi 
|[| where analytic expressions are derived from the minimal area law to first order in 
the inverse heavy quark mass. More recently, a relativistic QCD-string Lagrangian 
for spinless quarks has been deduced. 

Here we derive an effective Hamiltonian for a quark-antiquark system in the 
confinement region with the minimal area (MA) law of lattice gauge theory as input. 
In this arrangement, the gauge field's non-Abelian character is manifest through the 
minimal area law; complete non-Abelian expressions are therefore given only where 
useful, e.g., to clarify a point. Our method consists in expressing the usual canonical 
matrix elements from (|2]) in terms of Wilson loop expectation values derived from 
the MA law. In the effective theory these become first- quantized operators acting 
on state vectors. The results are fully relativistic. 
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2 Wilson loop expectation values 



We begin with the Lagrangian density assuming that the MA analysis has been 
carried out; hence fermion and gauge fields are interrelated at the outset, effectively 
reducing the degrees of freedom in (0). The idea is simple: gauge fields specified 
by MA relations are taken to be "external" to the Lagrangian; Euler-Lagrange non- 
linearities ([§|3]) are thereby hidden within undetermined matrix elements. 

From the energy-momentum tensor of (Q) the conserved quantities are the usual 

H = E / <$ x Qi h ■ H V " 9 K ) + 7 Vo + rnjjqj , (5) 

3 

J d 3 xq}(-iV)q 3 , (6) 

3 

J2 J d 3 xq}[x x (-iV) + , (7) 



P 
J 



3 

which with the MA law, 

i\n(W(C)) = aS mm , (8) 

define our bound state problem. Eventually, we take independent variations of (§) 
along the paths C. 

We consider the Wilson loop of a straight slice out of the qq minimal world sheet 
from time r', 

(W(C)) = ±(TrPexp(i9 <f dafA^x))) 

= i(TrPexp(-zs jH' dr [A ( Zl ) - z 1 ■ A( Zl ))) U(z' 2 , z[) 

2 dr[A {z 2 ) -z 2 - A(z 2 )})U(z lf ,z 2f )) . 

■2/ 

The U's here are straight-line path ordered exponentials, and the average is taken 
over gauge fields, 

JVAeM^[AMA] 
[ ' JVAexp(iS[A]) ' { ' 
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where the pure gauge action S[A] includes a gauge fixing term. 

The c-number path variables for a given gauge field (zj and Zj) are related to 
coordinate and mechanical-momentum matrix elements, 



maiZi = (qi(T;A)\ (p- gA(x)) \qi(T;A)) 
1 



Zj = Zj(r; A) = (q^r; A)\ x ^(r; A)) 



(10) 

(ii) 

(12) 



where is either a particle or antiparticle Euler-Lagrange solution (§-§) with a 
classical gauge field. These valence or quenched solutions prevent the possibility of 
time-backtracking not present in the MA picture. If x is the meson wavefunction 
in a product space of |gi) and | , a reasonable ansatz is the auxiliary condition 



where A± are positive or negative energy projection operators. In this approxima- 
tion, the Casimir operators reappear in the Hamiltonian ([5]). 

On the contrary, we point out that the quark degrees of freedom have not yet 
entered the dynamics except as a means to specify the physical gauge field configu- 
rations. The above identifications are made for later use. 

The enclosed area we parameterize in the Nambu-Goto form, 



with x^ = x m (t, a), x = t, x^ = -J^-, and x'^ = At the boundary, x(r, 1) = 
(zi(t)) and x(r, 0) = (z 2 (r)). In the usual straight-line and equal-times approxima- 
tions, which we assume here, the minimum is given by 



X - a« x aL 2) - a«xa: 




(13) 



s 




(14) 




a( Zl (r)) + (l-a)(z 2 (r)) . 



(15) 



5 



A general path variation of (|8]) for fixed endpoints and fixed time, 

J dr-(^ Tr P exp (tg j> dx^A^(x)) 



5(W(C)) -% 

"WW " WW)" 



r a a 
x g\ (5z 1 ■ — + 5zi • — ) [Aoizx) - zx ■ A(zi)} 
I ozi azi 

d d 

- (£z 2 -- h <Sz 2 • — ) A)(z 2 ) - z 2 ■ A(^ 2 ) 

oz 2 dz 2 



c?t / da 



dS . ... dS .... 
-Sx'" + -—dx 11 



dx'f"~ ' dx*"~ ' (H,) 

specified to 5z\ = <5z 2 , obtains in the center of momentum (see Appendix for details) 

((gA(z))) = — j— ^TrPexp^jfda^A^Cx))^^)^ 

/•<*> . . 

= a / dx-y ± x± ee p t , (17) 
Jo 

which are the desired Wilson loop expectation values with reference points chosen 
at the origin, A(0) = 0, and "_L" defined relative to the straight line, v_i_ ee v — 
(v ■ r)r. It should be clear that the gauge fields quantized in this average carry 
spatial dependence in the radial coordinate only. Evidently, this proceedure selects 
the physically realizable transverse polarizations. The time component is obtained 
by simple differentiation of (fflj with respect to r, 

±-i\n{W{C)) = —^—^I^TrPexp^^xM^x)) 

x g{A (z 1 ) - A (z 2 ) - [ Zl ■ A(*i) - z 2 • A(z 2 )]}^ 

^daS Xmin , (18) 



o 



yielding (using the proceedure described in the Appendix) 

<•<*> 



((gA (z))) = a [ * dxj ± = H t , (19) 
Jo 



and also 

■<*) 



((z x gA(z))) = a / dx(x x 7± x ± ) ee L t . (20) 
Jo 
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3 Relation to Flux Tube Model 



With the Wilson loop expectation values we re-express conserved quantities in 

terms of coordinate and velocity matrix elements. In the effective theory these are 

promoted to noncommuting quantum operators, observables requiring symmetriza- 

tion. We give the final expressions 
2 

i=l 
2 

i=l 
2 

J = [ X * X m ^i\ S yra + 

i=l 

This compares well with other results from the Wilson loop. Both spin and spin- 
independent Hamiltonians of for example, are reproduced on semi-relativistic 
reduction of (p2[). Equations ( ]21] - p3D in fact define the relativistic flux tube model 
as formulated by the present authors in 0, |3[. There, dynamical relations are derived 
from a heuristic Lagrangian in which the flux tube is fashioned as a simple constant 
energy density. The "tube operators" are equivalent to the Wilson loop expectation 
values above. This straightforward derivation serves to clarify the model's close 
relation to the underlying QCD. 



((gA Q {xi))) +Pm, 



(21) 



«<?A(x,))) 



sum 



(22) 



((x; x gA(xi))) 



u j 

+ — 

finrn. 



(23) 



4 Conclusion 

We have demonstrated that under the same assumptions as in previous work [III [71 
that the Wilson loop area law and the QCD Lagrangian yields a relativistically valid 
picture of dynamical confinement. We emphasize that we have not "solved" QCD, 
but that our ignorance can be distilled into three loop expectation values each of 
which is equal to that of the mechanical RFT model. We find that 

((gAo)) = H t , (24) 
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«0A» = Pt , 
((r x gA)) = L t , 



(25) 
(26) 



for the tube energy, momentum, and angular momentum respectively. In addition 
we have shown that the natural equation of motion is the Salpeter equation in which 
the "covariant tube substitution" || has been made. 

The verification of the RFT model structure promises more than to legitimize a 
physically reasonable model. The close relation of the RFT model to lattice QCD 
should shed light on both subjects. In addition, a systematic program of improving 
the RFT model to include field fluctuations || can now be envisioned. 



A Appendix 

Equation ( |TTD follows from the physical gauge field's spatial dependence in the 
straight line and equal times approximations 

^( r ) ~ exp(— ik ■ r) — > exp(— ikr) . (27) 
k k 

Then small angular variations of a given path at fixed time leave A unchanged. It 
will suffice to consider the transverse part of ([II]) with 

Szi(r;A) = 8 Zi (r;A') = 5z t , (28) 

for all gauge fields at a given r. Naturally, 

S(zi) = 5z { . (29) 

We write the Wilson loop in discrete form, 

1 N 

W{C) = -Tr P J] exp(igAt n [A (z ln ) - z ln ■ A(z ln )))U(z' 2 , z[) 

6 n=0 

x exp(-igAt n [Ao(z 2 n) - z 2 „ • K(z 2n )])U (z lf , z 2f ) , (30) 
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and take the path variations of (||) according to 



5{iln(W{C))) 



(W(C)) 



5(W(C)) 



N 



1 i 1 v Q 



(W{C))\i 







dzin 

d 



dz 



In 



+ bz 2n -j- + 5z 2n --^ )W(C)) 

OZ 2n OZ 2n I 



(31) 



giving equation flip]), 

1 r T f 1 / / 2 d d 

, —7—rr / dr—( TrPexp [%q <p dx^AJx)) y^(<5z; • — h 5z; ■ — — ) f(zi,Zi) 

(W{C))Jt> 3\ yyy Jc tk dz * dz i 



a l dr J da 



os d . u ds d . u 

Sx 1 " + -77— —SX^ 



dx 1 ^ da 



dx^ dr 



where 



-1)% ■ gA(zi) - gA (zi)) 



(32) 



(33) 



The temporal derivative is transferred from the variations by partial integration, so 
that 



(W(C)) Jr' ^13 

= a p dr Jo 1 da 



P drZ, |( TrPexp (zg § c dx»A,(x))(--f T £- + £) /(*, z 4 ) > ■ 5z 



,(34) 



yielding in the bz\ = 5z 2 (= 5zj_) case of interest 



E«4r/(^,*i)» 

i 2 



a Jo da (l^), 



(35) 



or 



II Kt \ Kf \\\ 1/ \ / \l + (1 - <t)(z 2 ±) 

((<,a( Zi ) - S A(, 2) » = «!<„> - ( z 2 )i y o ^ |1 _ (g(zii) + (1 _ g)(z2i))2|1/2 . 

This obtains (|ITD when reference points are chosen at the origin of coordinates 



(36) 
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Also, with the A reference point chosen at the origin, A (0) = 0, and fll8|) can 
be written 

1 1 TrP JdAtg[A t (z At )-z At -At(z At )]JVA^ t exp(ig§ c dx' i A IJ .)e X .p{iS[A]) 
(W{C)) 3 J VAexp(iS[A]) 

= a(z t )tida[l-a 2 (i t )i] 1 /\ (37) 

with 

JdA t z At JVA^ t exp(iS[A]) 

{Zt) ~ JVAeMiS{A}) • (38) 



Then, the functional derivative of fl37|) with respect to z^/ is 

1 | Tr P (-gAQ / PA^ exp (ig § c dx»A„) exp (ig[A]) 
(W(C)) 3 J 2X4exp(jS[A]) 

= afa> J.' ^)Wi[l - ^Wi]" /i J ^g|g ■ (39) 
Taking the scalar product of the above with z^, and integrating over A' t , gives 

{(z-gA)) = a f" dx±\ 1± , (40) 
Jo 

which is used in obtaining fll9"|). Taking the vector product of ( |39"D with z^j yields 
after integration over A' t , 

/■<*) 

((z x gfAf^))) = a dxCx x 7_lx_l) • (41) 
Jo 
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